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Proper  water  management  (sufficient  membrane  hydration  and  effective  liquid-water  removal  from  cat¬ 
alyst  layers,  gas-diffusion  layers,  and  channels)  is  key  to  achieving  and  maintaining  the  high  performance 
of  polymer  electrolyte  fuel  cells  (PEFCs).  Specifically,  liquid  droplet  formation  and  removal  is  an  impor¬ 
tant  issue  in  water  management  for  PEFCs.  In  the  first  part  of  this  series,  explicit  analytical  solutions 
were  derived  by  approximating  the  droplet  shape  as  spherical  and  the  Navier-Stokes  equations  were 
numerically  solved  to  compute  the  forces  over  a  droplet.  In  this  second  part  of  the  series,  results  from 
numerical  simulation  using  the  Volume-of-Fluid  (VOF)  method  and  experiments  are  reported  and  they 
are  compared  with  the  analytical  solutions  derived  previously.  More  specifically,  experimental  visual¬ 
ization  of  droplet  dynamics  in  a  micro  channel  is  carried  out.  A  high  resolution  CCD  camera  is  employed 
to  capture  the  droplet  shape-change  and  detachment.  Extensive  numerical  simulations  via  VOF  are  also 
performed  to  investigate  droplet  dynamics.  Reasonably  good  agreements  between  analytical  solution 
and  VOF  numerical  simulation  and  experiment  are  obtained. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Proper  water  management  (sufficient  membrane  hydration  and 
effective  liquid-water  removal  from  catalyst  layers,  gas-diffusion 
layers,  and  channels)  is  key  to  achieving  and  maintaining  the  high 
performance  of  polymer  electrolyte  fuel  cells  (PEFCs).  Specifically, 
liquid  water  removal  from  PEFC  channels  is  an  active  and  impor¬ 
tant  area  of  research  in  PEFCs.  Specifically,  water  droplet  formation 
frequently  occurs  at  the  gas  diffusion  layer  (GDL)  surfaces  in  the 
gas  flow  channel  of  PEFCs  [1-3].  One  source  of  liquid  water  is 
from  condensation:  when  the  vapor  partial  pressure  reaches  the 
saturated  one,  excessive  water  in  the  vapor  phase  condenses  to 
liquid  water.  The  other  is  from  liquid-water  production  directly 
by  the  electrochemical  reaction  in  catalyst  layers.  Liquid  water 
produced  in  the  cathode  catalyst  layer  can  be  absorbed  into  mem¬ 
brane  and  transported  to  anode  side  by  diffusion  or  permeation. 
It  can  also  be  delivered  to  the  cathode  gas  channel  via  the  GDL. 
At  the  hydrophobic  GDL  interface,  liquid  can  aggregate  to  form 
droplets  and  grow  its  volume  until  it  is  removed  by  the  reactant 
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gas  flow.  Droplet  formation  can  increase  transport  polarization, 
and  consequently  reduce  cell  performance. 

Droplet  dynamics  in  a  micro/mini  channel  have  been  studied 
through  simplified  analysis,  numerical  simulation,  and  visualiza¬ 
tion  experiments.  Early  studies  on  droplet  dynamics  for  PFEC 
application  mainly  used  optical  visualization  to  examine  droplet 
formation  and  removal.  Yang  et  al.  [4]  placed  a  transparent  win¬ 
dow  on  a  multi-channel  PEFC,  showing  droplet  formation  and 
channel  blocking.  Wang  et  al.  [5]  showed  experimental  visualiza¬ 
tion  of  two-phase  flow  in  multiple  gas  flow  channels,  indicating 
droplet  formation.  They  also  linked  the  two-phase  flow  to  the 
pressure  drop  and  further  cell  voltage  oscillation  when  decreasing 
the  stoichiometry  to  around  2.0.  Chen  et  al.  [6]  conducted  experi¬ 
ments  on  an  operating  PEFC  with  a  transparent  channel  wall  and 
was  able  to  observe  the  formation  and  breakup/detachment  of 
water  droplets  on  the  GDL  surface.  They  also  pioneered  droplet 
instability/detachment  visualization  experiments  using  a  specially 
designed  visualization  apparatus  with  transparent  channel  walls, 
and  developed  droplet-instability  windows.  Ma  et  al.  [7]  linked 
the  droplets  in  channels  to  the  pressure  drop  and  use  it  for  water 
removal  diagnosis.  Wang  [8]  indicated  that  droplets  form  and 
attach  to  the  side  wall  of  a  channel.  Theodorakakos  et  al.  [9] 
used  transparent  PEFCs  to  study  droplet  detachment.  Static  and 
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Nomenclature 

A  area(m2) 

D/d  droplet  diameter  (m,  mm) 

Dh  hydraulic  diameter  (m) 

F  force  (N) 

H  channel  height  (m) 

i  direction  vector  ( - ) 

h  normal  vector  (-) 

p  pressure  (pa) 

R  local  curvature  (m) 

RH  relative  humidity  (-) 

Sm  mass  source  term  (kg  m-3  s-1 ) 

r  droplet  radius  (m) 

t  tangential  vector  ( - ) 

u  velocity  (m  s-1 ) 

W  channel  width  (m) 

Wer  Weber  number  (-) 

Greek  letters 

a  volume  fraction  (-) 

y  surface  tension  (Nm-1) 

0  contact  angle  (rad,  deg) 

/a  viscosity  (kg  nrr 1  s-1 ) 

p  density  (kg nrr3) 


dynamic  contact  angles  were  measured  and  compared  with  VOF 
(Volume  of  Fluid)  simulation  results.  The  VOF  has  been  regarded 
as  a  useful  numerical  tool  to  study  the  droplet  dynamics  thanks 
to  its  capability  of  tracking  liquid-gas  interface.  Zhu  et  al.  [10] 
used  VOF  for  parametric  study  on  the  impact  of  channel  dimen¬ 
sion  on  droplet  detachment  velocity.  Bazylak  et  al.  [11]  studied 
droplet  emergence  and  detachment  through  VOF  simulation  and 
ex  situ  experiment.  Hao  and  Cheng  [12]  analytically  predicted 
dynamic  contact  angle  and  droplet  velocity  in  a  micro  channel. 
The  results  are  compared  to  LBM  (Lattice  Boltzmann  Method) 
simulation.  An  explicit  analytical  model  to  predict  the  droplet  insta¬ 
bility  in  the  viscous-force-dominating  flow  regime  (i.e.,  Reynolds 
number  is  vanishingly  small)  was  proposed  by  Chen  et  al.  [6] 
using  the  force  balance  approach.  As  mentioned  earlier,  visual¬ 
ization  experiments  were  also  conducted  to  generate  data  for 
comparison  with  model  prediction.  Kumbur  et  al.  [13]  proposed 
similar  approach  along  with  experimental  work.  More  recently, 
Chen  [18]  presented  an  explicit  analytical  model  for  predict¬ 
ing  droplet-detachment  velocity  in  the  inertia-dominating  flow 
regime. 

Though  many  studies  have  been  proposed  to  understand  droplet 
dynamics,  explicit  analytical  models  that  can  be  employed  as 
tools  for  PEFC  design  and  optimization  are  still  in  great  need 
and  their  comparisons  with  experimental  data  are  lacking.  In  the 
first  part  of  this  paper  series,  analytical  solutions  for  predicting 
droplet  deformation  and  detachment  were  derived  by  approxi¬ 
mating  the  droplet  shape  as  spherical.  The  Weber  number  and 
Reynolds  number  are  linked  to  droplet  shape-change  and  removal 
through  analytical  solutions.  In  the  present  work  (which  is  the  sec¬ 
ond  part),  experiments  and  VOF  numerical  simulation  (without 
approximating  the  droplet  shape  as  spherical)  on  droplet  defor¬ 
mation  and  detachment  were  conducted  and  the  results  were 
compared  with  the  analytical  solutions.  We  hope  the  discussion 
of  forces,  deformation,  and  detachment  will  provide  a  compre¬ 
hensive  picture  of  droplet  dynamics  over  GDL  surfaces.  Again 
the  results  are  beneficial  to  the  fundamental  understanding  of 
droplet  formation  and  detachment  at  GDL  surface,  which  is  of 
vital  importance  to  flow  field  design  and  water  management  for 


Fig.  1.  Schematic  and  picture  of  experimental  setting  (top);  and  schematic  of  a 
droplet  in  a  gas  flow  channel  (bottom). 


PEFCs.  For  example,  the  analytical  solutions  contain  parameters 
such  as  the  channel  dimension,  surface  tension,  and  detachment 
velocity,  therefore  they  can  be  used  to  optimize  flow  channels  and 
droplet  removal  strategy.  Further,  droplet  dynamics  over  GDL  sur¬ 
faces  is  one  major  area  of  research  in  channel  two-phase  flow. 
There  are  other  areas  such  as  liquid  flow  over  the  hydrophilic  chan¬ 
nel  sidewall  and  droplet  removal  when  it  reaches  the  sidewalls. 
It  is  also  important  to  investigate  all  the  issues  by  consider¬ 
ing  the  entire  channel  due  to  the  inherent  connection  between 
them.  This  will  allow  the  material  nature  of  the  flow  field  as 
a  whole  to  be  accounted  for  in  studying  channel  two-phase 
flow. 


2.  Analytical  solutions 


In  the  first  part  [14],  the  drag  coefficient  of  a  droplet  in  a  micro 
channel  as  a  function  of  the  Reynolds  number  is  generalized  using: 

CD  =  axRebH  (1) 

where  a  =  46.247  x  (D/H)01 757  and  b  =  0.21 58  x  (D/H)  -  0.6384  in 
the  fully  developed  region  for  the  channel  considered.  With  the 
above  generalized  formula,  the  detachment  velocity  in  terms  of  the 
Weber  number  ( Wer  =  pu2R0ly)  can  be  expressed  by: 


Wer 


4  7 r  sin2  0S  sin(l /2 ){0a  -  0r) 
a  ( 0S  -  sin  0S  cos  0S) 


Re-f 


(2) 


Droplet  deformation  due  to  its  obstruction  can  be  evaluated 
by  the  Young-Laplace  equation  and  Bernoulli’s  equation,  which  is 
quantified  by  the  change  of  local  curvature  radius: 


dR 

~R 


f  1  -  ^(0S  -  COS  0S  sin  Os) 


(3) 


where  R,  R0  and  Ac  are  radius  of  local  curvature,  initial  radius  of 
droplet  and  cross  sectional  area  of  channel,  respectively.  Details 
of  the  analytical  model  can  be  found  in  the  first  part  of  this  paper 
series  [14]. 


Table  1 

List  of  experiment  cases. 
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Case  No. 

GDL 

Channel  cross  section 
dimension  (mm2) 

Contact  angle 

f^drop/f^H 

Air  flow  velocity  (m  s-1 ) 

Flow  condition 

1 

Carbon  paper  (wet  proof  30%) 

1.6  x  1.0 

128° 

0.90 

0.0-5.21 

Fully  developed 

2 

PTFE  tape 

1.6  x  1.0 

115° 

0.89 

0.0-5.21 

Fully  developed 

3 

Carbon  paper  (wet  proof  30%) 

1.6  x  1.0 

128° 

0.69-1.01 

0.0-8.75 

Entrance  length 

4 

Carbon  paper  (wet  proof  30%) 

1.6  x  1.0 

128° 

0.69-0.98 

0.0-8.75 

Fully  developed 

3.  Experimental 

Fig.  1  shows  the  experimental  setup  that  was  specially  designed, 
fabricated  and  employed  to  study  droplet  dynamics  in  micro/mini 
channels.  The  gas  flow  channels  were  fabricated  within  both  Alu¬ 
minum  and  plastic  thin  plates,  and  both  were  used  in  experiment 
to  ensure  independence  of  the  experimental  data  on  channel 
wall  materials.  The  channel  has  a  cross-sectional  dimension  of 
1.0  mm  x  1.6  mm  with  length  of  4.0  cm.  The  GDL  material  is  the 
Toray  carbon  paper  with  30%  PTFE  loading.  A  transparent  plate 
was  placed  on  the  channel  plate  to  provide  optical  window.  A  CCD 
camera  was  mounted  over  a  microscope  for  visualizing  droplet 
dynamics.  The  base  plate  was  fabricated  by  CNC  machining  on  Alu¬ 
minum  plates  where  air  and  water  injection  fittings  were  mounted. 
The  hydrophobic  carbon  paper  is  placed  on  the  base  plate  and 
clamped  by  the  channel  and  the  top  transparent  plates.  Liquid 
droplet  is  generated  by  supplying  water  from  the  base  plate  via  a 
small  hole  through  the  carbon  paper.  Water  permeates  the  porous 
layer  and  forms  a  droplet  on  the  GDL  surface.  Since  the  contact  angle 
of  hydrophobic  materials  are  104-128°,  the  droplet  diameter  can 
reach  over  channel  height  ( 1  mm)  without  touching  the  top  surface. 
Once  a  droplet  is  formed  in  a  static  flow,  the  diameter  of  droplet 
is  measured  and  gas  velocity  is  increased  gradually  by  0.1  ms-1 
until  the  droplet  is  removed  from  the  GDL  surface.  In  real  fuel  cell 
operation,  however,  liquid  water  is  added  to  form  a  droplet  in  air 
flow  environment  till  its  detachment.  This  likely  makes  no  differ¬ 
ence  in  terms  of  the  forces,  deformation  and  detachment  of  droplets 
because  it  is  apparent  that  the  physical  presence  of  liquid  droplet  is 
the  major  cause.  Further,  the  present  method  allows  us  to  precisely 
measure  the  droplet  size  because  it  remains  perfect  spherical  shape. 
When  there  is  air  flow,  the  droplet  shape  can  change,  making  it  dif¬ 
ficult  to  measure  the  droplet  dimension.  As  stated  before,  the  water 
of  droplets  may  come  from  vapor  condensation  in  supersaturated 
gas  or  liquid  transport  through  GDLs.  Because  of  low  temperature, 
condensation  occurring  at  droplet  surface  is  negligible,  therefore 
it  is  the  latter  source  that  is  primarily  considered  in  the  experi¬ 
ment.  However,  the  sources  of  water  have  little  effect  on  the  forces, 
deformation,  and  detachment  of  droplets  because  they  are  primar¬ 
ily  determined  by  the  physical  presence  of  liquid  droplets.  Similar 


strategy  is  also  used  in  Ref.  [6].  The  water  injection  hole  was  drilled 
on  the  base  plate  about  3.75  cm  away  from  the  channel  manifold,  i.e. 
the  flow  near  the  droplet  site  is  in  the  fully  developed  condition.  By 
placing  honeycomb  (porous  material)  in  the  upstream  gas  channel 
near  the  droplet  inlet  hole,  the  flow  field  at  the  droplet  site  can  be 
modified,  i.e.  the  droplet  is  subject  to  air  flow  in  the  entrance  region. 
The  entire  test  section  was  placed  on  the  pneumatic  isolated  table 
to  eliminate  any  noises  from  external  disturbance.  Air  flow  rate  is 
controlled  by  the  mass  flow  controller  in  a  fuel  cell  test  station  with 
100%  relative  humidity.  Details  of  the  experimental  cases  are  listed 
in  Table  1.  Droplets  on  the  surfaces  of  the  two  materials  are  shown 
in  Fig.  2.  Note  that  the  experimental  set  up  is  similar  to  that  used  by 
Adroher  and  Wang  [15],  which  was  designed  to  study  two-phase 
flow  pattern  and  pressure  drop. 

4.  Numerical  simulation:  Volume  of  Fluid  ( VOF)  method 

A  three-dimensional  model,  based  on  the  VOF  method  [16] 
(which  is  one  of  the  most  popular  numerical  techniques  to  trace 
the  interface  of  two  phases  and  study  multiphase  flows),  was 
developed  and  used  to  simulate  the  water-droplet  dynamics  (defor¬ 
mation,  detachment).  With  this  method,  the  equations  that  govern 
mass  and  momentum  conservations  are  solved  in  conjunction  with 
the  volume  fraction  equation.  The  former  can  be  expressed  as: 

9 

7^p  +  V-(pu)  =  Sm  (4) 

7^(pu)  +  V  •  (put)  =  - Vp  +  V  •  (/x(Vu  +  Vu7))  +  pg  +  FVol  (5) 

The  momentum  source  term  FVoi  on  the  right  hand  side  of  the 
equation  takes  into  account  the  surface  tension  y,  surface  curvature 
kk  and  gradient  of  the  volume  fraction  Vak.  The  volume  fraction 
ak  is  defined  as  the  portion  of  a  computational  cell  filled  with  kth 
phase  of  fluid.  In  this  study,  the  primary  fluid  is  set  as  air  and  the 
secondary  fluid  is  liquid  water.  Thus,  any  computational  cell  filled 
with  only  air  has  the  value  of  zero  for  awater  •  For  the  computational 
cell  filled  with  liquid  water,  awater  becomes  unity.  A  value  between 
zero  and  unity  means  the  cell  contains  the  two-phase  interface.  To 


Fig.  2.  The  contact  angles  of  two  hydrophobic  materials. 
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satisfy  the  continuum  theory  of  fluid,  the  summation  of  a  in  a  cell 
satisfies: 


Oi  water  +  Ofair  —  1 


(6) 


The  volume  fraction  is  governed  by  the  following  volume  frac¬ 
tion  equation: 

9 

air  Pair)  +  ^  '  i^air  Pair  ^  air)  —  0  (7) 


This  volume  fraction  conservation  equation  is  correlated  with 
the  mass  and  momentum  equations  through  the  volumetric  force 
term  (FVo/)  and  solved  every  iteration  step  for  all  computational 
cells.  The  properties  of  fluids  at  the  interface  between  gas  and  liquid 
are  weight-averaged  by  the  volume  fraction  in  the  cell,  for  example: 


P  —  Pair  +  &  water  (Pwater  ~  Pair) 


(8) 


Through  this  averaging,  fluids  in  interface  can  be  regarded  as 
single  phase  and  only  one  set  of  governing  equations  needs  to  be 
solved  in  the  interface  region.  For  the  local  curvature  of  interface, 
the  Young-Laplace  equation  can  be  utilized  considering  two  radii 
in  the  orthogonal  direction  ( R !  &  R2 ): 

Ap  =  y{r^  +  rT)  (9) 

The  surface  tension  force  is  related  to  the  momentum  source  by 
using  the  CSF  (continuum  surface  force)  model  [17],  in  which  the 
force  can  be  expressed  as: 


Fvoi  =  Y 


pkkVak 
( V2)(Pfe  +  Pj ) 


(10) 


where,  kk  =  V  •  n,  n  =  V  ak,  and  h  =  n/|n|. 

At  the  contact  line  on  the  GDL  surface,  droplet  and  air,  the  surface 
normal  h  is  determined  by  the  static  contact  angle  0S : 


n  =  nw  cos  0S  +  tw  sin  0S 


(11) 


where  hw  and  tw  are  the  normal  and  tangential  vectors,  respec¬ 
tively,  to  the  wall  surface.  Based  on  this  surface  normal,  the  local 
curvature  near  the  wall  boundary  and  surface  tension  force  are 
determined  in  turn.  In  the  present  work,  about  576,000  hexahe- 
dral  computational  cells  were  used  in  the  VOF  simulation.  For  the 
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Fig.  4.  Air  pressure  distribution  near  a  water  droplet  (top)  and  droplet  shape  change 
(bottom). 


S.C.  Cho  et  al. /Journal  of  Power  Sources  210  (2012)  191-197 


195 


Fig.  5.  Comparison  of  analytical  solution  with  the  VOF  results. 


pressure  &  velocity  coupling,  the  SIMPLE  method  is  used,  and  the 
governing  equations  are  discretized  through  PRESTO!  using  the 
second  order  upwind  method.  In  this  work,  commercial  CFD  (Com¬ 
putational  Fluid  Dynamics)  solver  Fluent  (version  6.3.26)  was  used 
to  carry  out  the  VOF  numerical  simulation. 

5.  Results  and  discussion 

5.1.  Droplet  shape  change 

Fig.  3  shows  an  example  of  the  water  droplet  dynamics  sim¬ 
ulated  using  the  3D  VOF  method.  A  droplet  forms  at  the  site  of 
the  hole  by  liquid  injection.  At  low  velocity,  the  wall  adhesion  is 
sufficient  to  keep  the  droplet  still  on  the  surface.  As  the  gas  flow 
rate  increases,  the  drag  force  increases  rapidly  [14],  which  eventu¬ 
ally  overcomes  the  adhesion  force  and  removes  the  droplet  from 
the  surface.  In  the  first  part  of  this  paper  series  [14],  the  droplet 
shape  change  caused  by  the  pressure  difference  is  quantified  by  an 
analytical  solution  derived  by  approximating  the  droplet  shape  as 
spherical.  Fig.  4  shows  the  gas  pressure  contours  predicted  by  the 
VOF  simulation,  indicating  the  gas  phase  pressure  can  vary  substan¬ 
tially  due  to  the  obstruction  of  the  droplet.  At  the  stagnation  point 
in  the  front  of  the  droplet,  the  local  gas  pressure  is  increased,  while 


Fig.  6.  Visualization  of  droplet  dynamics  in  a  micro  gas  flow  channel. 


at  the  location  where  the  cross  sectional  area  for  gas  flow  decreases 
to  the  minimum  and  the  local  gas  velocity  increases  to  maximum, 
local  pressure  becomes  low.  It  can  be  seen  that  the  droplet  shape 
is  changed  and  slightly  deviates  from  the  spherical  shape:  at  the 
lower  gas  pressure  site  the  surface  curvature  is  larger  to  balance 
the  pressure  difference  between  gas  and  liquid;  at  the  higher  gas 
pressure  site  the  liquid  withdraws  resulting  in  a  smaller  curvature. 

The  droplet  deformation  caused  by  the  pressure  variation  can 
be  directly  calculated  through  the  liquid-gas  interface  shape  pre¬ 
dicted  by  the  VOF  method,  and  the  result  from  a  case  study  is  plotted 
in  Fig.  5  together  with  that  predicted  by  the  analytical  solution.  For 
both  2D  and  3D  simulation,  the  VOF  results  show  slightly  larger 


Fig.  7.  Experimental  observation  of  droplet  shape  change  on  hydrophobic  GDL  surface,  PTFE  (case  2)  (0S  =  115°). 


196 


5.C.  Cho  et  al. /Journal  of  Power  Sources  210  (2012)  191-197 


Fig.  8.  Comparison  of  droplet  deformation  between  the  analytical  solution  and 
experimental  data. 

deformation  but  the  trends  are  in  good  agreement.  The  discrepancy 
may  be  due  to  that  the  average  velocity  is  used  in  the  analytical  solu¬ 
tion  whereas  there  is  a  local  velocity  variation  in  the  VOF  numerical 
simulation. 

Experimentally,  the  droplet  shape-change  was  captured  using  a 
high-resolution  camera.  GDL  surfaces  were  examined  in  this  work, 
the  carbon  paper  (0  =  128°)  and  PTFE  surface  (0  =  104°),  see  Fig.  2. 
Fig.  6  displays  the  visualization  of  droplet  emerging,  growth,  and 
detachment  on  a  carbon  paper  surface  whereas  Fig.  7  shows  the 
visualization  of  droplet  deformation  on  the  PTFE  surface.  In  Fig.  6, 
the  droplet  was  introduced  into  the  channel  where  exists  air  flow, 
similar  to  what  happens  in  PEFC  channels.  Fig.  8  compares  the 
experimental  data  with  the  analytical  solution.  Due  to  the  hystere¬ 
sis  of  the  contact  angle,  it  was  difficult  to  determine  the  droplet 
front  line  (or  upstream  front)  optically  at  high  velocity  in  experi¬ 
ment;  consequently,  we  focused  on  the  regime  with  velocity  less 
than  2.0  m  s-1 .  It  can  be  seen  that  acceptable  agreement  is  achieved 
-  the  observed  deformation  is  less  than  5%  for  all  cases  consid¬ 
ered  in  the  experiment.  Thus,  droplet  deformation  can  be  negligible 
in  the  velocity  regime  considered  in  the  present  work,  and  its 
effect  on  the  flow  drag  coefficient  can  be  neglected.  Large  defor¬ 
mation  is  observed  experimentally  at  high  air  velocity,  however,  it 
is  extremely  challenging  to  quantify  this  due  to  the  hysteresis  of 
contact  angle  as  mentioned  previously. 

5.2.  Detachment  velocity 

The  analytical  results  of  the  detachment  velocity  using  the  new 
correlation  for  the  drag  coefficient  (Eqs.  (1)  and  (2))  are  compared 
to  various  experimental  data  as  presented  in  Fig.  9.  Because  the  drag 
coefficient  in  the  entrance  region  varies  with  the  distance  from  the 
inlet,  we  skip  the  analytical  solution  for  that  in  the  entrance  region. 
In  Fig.  9,  the  experimental  data  is  compared  with  analytical  solution 
for  1 .0  mm  x  1 .6  mm  channel  with  hydrophobic  carbon  paper  (wet 
proofed  with  30  wt%  PTFE  loading).  The  contact-angle  hysteresis  is 
set  as  40°  in  the  analytical  solution.  The  operating  condition  is  1  atm 
and  25  °C.  It  shows  a  good  agreement  with  the  new  correlation  of 
the  drag  coefficient  for  the  fully  developed  region.  In  the  experi¬ 
ment,  the  detachment  velocity  for  droplets  in  the  entrance  length 
is  higher  than  that  in  fully  developed  condition,  which  is  consis¬ 
tent  with  the  force-balance  analysis  in  [14].  Note  that  the  droplets 
considered  experimentally  are  relatively  large  with  diameter  over 
0.5  mm. 


From  the  analytical  solution,  a  Re-We  correlation  was  devel¬ 
oped  at  the  detachment  velocity  [14],  which  is  plotted  in  Fig.  10  in 
comparison  with  the  experimental  data.  The  Re-We  relation  based 
on  new  correlation  of  the  drag  coefficient  (i.e.  Eq.  (2))  is  shown  for 
various  H/r  values.  For  the  comparison,  H/r  values  of  several  points 
are  indicated  in  the  figure.  The  ratio  between  droplet  radius  and 
channel  height  (H/r)  of  experimental  data  ranges  from  1.9  to  2.9. 
The  experimental  Weber  numbers  show  good  agreement  with  the 
analytical  solution. 

The  analytical  solution  is  also  compared  to  other  experimental 
work  in  larger  channel  dimension  (7.0  mm  x  2.7  mm)  with  similar 
GDL  materials  (hydrophobic  carbon  paper  &  cloth)  [9].  The  aver¬ 
age  static  contact  angle  in  the  experiment  is  130°  for  carbon  paper 
and  145°  for  carbon  cloth  [9,18].  In  Fig.  11,  the  experimental  data 
is  compared  with  analytical  results.  The  present  analytical  model 
shows  acceptable  agreement  with  experimental  data.  In  this  study, 
the  carbon  paper  has  a  higher  detachment  velocity  than  the  car¬ 
bon  cloth  and  the  analytical  model  also  follows  the  same  trend.  As 
expected,  the  detachment  velocity  decreases  as  droplet  diameter 
increases. 


Fig.  10.  Comparison  of  the  Weber  number  between  the  experiment  data  and  ana¬ 
lytic  solutions. 
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Fig.  11.  Comparison  of  the  detachment  velocity  between  analytic  solutions  and 
experimental  data  [9]. 


Further,  in  PEFCs  various  reactants  and  flow  conditions 
are  encountered,  such  as  different  temperatures  and  reactants 
(hydrogen,  air,  and  reformate).  For  example,  PEFCs  usually  work  in  a 
range  of  temperature  such  as  between  20  and  80  °C.  The  surface  ten¬ 
sion  will  change  from  around  0.0728  N  m-1  at  20  °C  to  0.0626  N  m-1 
at  80  °C,  about  15%  difference.  Gas  composition  changes  with  tem¬ 
perature  as  well.  The  analytical  solution,  however,  can  apply  to 
those  conditions  and  different  gases  by  considering  the  correspond¬ 
ing  properties  such  as  surface  tension  and  gas  density.  For  example, 
the  change  in  the  detachment  velocity  between  20  and  80  °C  can 
be  directly  evaluated  using  Eq.  (2). 

Again,  understanding  the  detachment  velocity  as  a  function  of 
material  properties,  channel  dimension,  and  operating  condition  is 
important  to  optimize  the  GDL  surface  properties,  water  manage¬ 
ment,  and  control  strategy.  Droplets  on  GDL  surfaces  cause  liquid 
accumulation  and  can  block  channel  flows  when  several  droplets 
merge.  By  choosing  proper  material  properties  (e.g.  contact  angle, 
see  the  Weber  number)  and  channel  dimension  using  the  analyt¬ 
ical  solution,  the  detachment  velocity  can  be  reduced  and  droplet 
effect  can  be  diminished. 


6.  Conclusion 

In  the  present  work,  the  analytical  solutions  (obtained  by 
approximating  the  water  droplet  shape  as  spherical)  for  droplet 
deformation  and  detachment  in  micro  gas  flow  channels  were  com¬ 
pared  with  the  VOF  numerical  simulation  and  experimental  data. 
We  used  both  the  VOF  and  experimental  visualization  to  evaluate 
the  dynamics  of  droplet  deformation.  The  pressure  distributions 
around  liquid  water  droplet  computed  numerically  indicates  a  high 
pressure  near  the  stagnation  point  and  low  pressure  at  the  site 
where  the  air  flow  accelerates.  The  pressure  difference  changes  the 
curvature  of  droplets,  but  the  droplet  deformation  is  small  at  rela¬ 
tively  low  velocities  (<2  m  s-1 ).  Quantitative  comparisons  between 
the  analytical  solution  and  experiment/numerical  results  show 
reasonably  good  agreement.  Detachment  velocity  measured  exper¬ 
imentally  is  also  compared  with  the  analytical  solution  -  acceptable 
agreement  with  experimental  data  for  various  channel  geometries 
is  achieved.  Lastly,  the  We-Re  correlation  at  detachment  veloc¬ 
ity  derived  from  the  analytical  solution  was  also  compared  to  the 
experimental  data,  indicative  of  good  agreement  in  the  range  of  H/r 
values  considered  in  the  present  work.  The  fundamental  knowledge 
obtained  in  the  present  work  is  important  to  the  understanding  of 
the  droplet  dynamics  in  the  gas  flow  channels  of  PEFCs  and  the 
optimization  of  GDL/channel  design  and  operating  conditions. 
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